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Abstract 

|We study the operator theory associated with such infinite graphs G 
as occur in electrical networks, in fractals, in statistical mechanics, and 
even in internet search engines. Our emphasis is on the determination 
of spectral data for a natural Laplace operator associated with the 
graph in question. This operator A will depend not only on G, but 
also on a prescribed positive real valued function c defined on the edges 
in G. In electrical network models, this function c will determine a 
conductance number for each edge. We show that the corresponding 
Laplace operator A is automatically essential selfadjoint. By this we 
mean that A is defined on the dense subspace V (of all the real valued 
functions on the set of vertices G° with finite support) in the Hubert 
space l 2 (G°). The conclusion is that the closure of the operator A is 
selfadjoint in Z 2 (G°), and so in particular that it has a unique spectral 
resolution, determined by a projection valued measure on the Borel 
subsets of the infinite half-line. We prove that generically our graph 
Laplace operator A = A c will have continuous spectrum. For a given 
infinite graph G with conductance function c, we set up a system of 
finite graphs with periodic boundary conditions such the finite spectra, 
for an ascending family of finite graphs, will have the Laplace operator 
for G as its limit. 
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1 Introduction 



The infinite graphs we consider live on a fixed countable infinite set, say L. 
Starting with such a set L (subject to certain axioms, listed below), we get 
a notion of edges as follows: Select distinguished pairs of points in L, say x 
and y, and connect them by a "line," called edge. In physics, when a vertex 
x is given, the set of vertices connected to x with one "edge" is called a set 
of neighbors, or nearest neighbors. Initially we do not assign direction to 
the edges. So, as it stands, an edge e is defined as a special subset {x, y} 
for selected points x, y in L. Think "nearest" neighbors! 

A direction is only assigned when we also introduce a function I on 
edges e, and then this function / is assumed to satisfy I(x,y) = —I(y,x). 
In electrical networks, such a function / may represent a current induced 
by a potential which is introduced on a graph with fixed resistors. So only 
if a current function / is introduced can we define a direction to edges, 
as follows: We specify source s(e) = x, and terminal vertex t(e) = y if 
I(x, y) > 0. meaning that the current flows from x to y. 

In this paper we study the operator theory of infinite graphs G, with 
special emphasis on a natural Laplace operator associated with the graph in 
question. This operator will depend not only on G, but also on a positive real 
valued function c defined on the edges in G. In electrical network models, 
the function c will determine a conductance number for each edge e. If 
e = (xy) connects vertices x and y in G, the number c(e) is the reciprocal of 
the resistance between x and y. Hence prescribing a conductance leads to 
classes of admissible flows in G. When they are determined from Ohm's law, 
and the Kirchhoff laws, it leads to a measure of energy, and to an associated 
graph Laplacian. We identify the Hilbert space 7i{G) which offers a useful 
spectral theory, and our main result is a theorem to the effect that the 
graph Laplacian is essentially selfadjoint, i.e., that its operator closure is a 
selfadjoint operator in 7i{G). 

Let G = (G°, G 1 ) be an infinite graph, G° for vertices, and G 1 for edges. 
Every x in G° is connected to a set nbh(x) of other vertices by a finite 
number of edges, but points in nbh(x) are different from x; i.e., we assume 
that x itself is excluded from nbh(x); i.e., no x in G° can be connected to 
itself with a single edge. Let c be a conductance function defined on G 1 . 

Initially, the graph G will not be directed, but when a conductance is 
fixed, and we study induced current flows, then these current flows will give 
a direction to the edges in G. But the edges in G itself do not come with 
an intrinsic direction. 

We show that the Laplace operator A = A c is automatically essentially 
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selfadjoint. By this we mean that A is defined on the dense subspace T> 
(of all the real valued functions on G° with finite support) in the Hilbert 
space H = TL{G):= l 2 (G°). The explicit formula for the graph Laplacian 
A = A(g jC ) is given in (|3.6p in section [3] below which also discusses the 
appropriate Hilbert spaces. The conclusion is that the closure of the operator 
A is selfadjoint in Ti., and so in particular that it has a unique spectral 
resolution, determined by a projection valued measure on the Borel subsets 
of the infinite half-line R+; i.e., the spectral measure takes values in the 
projections in the Hilbert space: = l 2 (G°). We work out the measure. 

In contrast, we note that the corresponding Laplace operator in the con- 
tinuous case is not essentially selfadjoint. This can be illustrated for example 
with A = — (d/dx) 2 on the domain T> of consisting of all C 2 -functions on the 
infinite half-line IR+ which vanish with their derivatives at the end points. 
Then the Hilbert space is L 2 (1R + ). 

So our graph theorem is an instance where the analogy between the 
continuous case and the discrete case breaks down. 

A second intrinsic issue for the operator theory of infinite graphs G, is 
that generically our graph Laplace operator A = A c will have continuous 
spectrum. We prove this by identifying a covariance system which implies 
that the spectrum of the corresponding Laplace operator will in fact be 
absolute continuous with respect to Lebesgue measure on the half-line. 

In a third theorem, for a given infinite graph G with conductance func- 
tion c, we set up a system of finite graphs with periodic boundary conditions 
such the finite spectra, for an ascending family of finite graphs, will have 
the Laplace operator for G as its limit. 



2 Assumptions 

In order to do computations and potential theory on infinite graphs G, it 
has been useful to generalize the continuous Laplacian A from Riemannian 
geometry |AC04j to a discrete setting [BHS05] . [CS07j, |Kig03| , [HKK02] . 



However the infinities for graphs suggest an analogy to non-compact Rie- 
mannian manifolds, or manifolds with boundary. 

Once the graph Laplacian is made precise as a selfadjoint operator it 
makes sense to ask for exact formulas for the spectrum of A. Our Laplace 
operator A = A c is associated with a fixed system (G, c) where vertices and 
edges are specified as usual, G = (G^°\ G^); and with a fixed conductance 
function c : G^ 1 ' — > R+. See (I3.6P below for a formula. 

And as usual our Laplace operator, A = A c is densely defined in the 



4 



Hilbert space £ 2 (G^) of all square-summable sequences on the vertices of 
G; and if G is infinite, A c is not denned everywhere in £ 2 , but rather it has 
a dense domain T> in £ 2 . We show in the next section that A c is essentially 
selfadjoint for all choices of conductance function c. 

By a graph G we mean a set G^ of vertices, and a set G^ of edges. 
Edges e consist of pairs s,y G G^°\ We write e = (xy); and if (xy) G G^ 
we say that x ~ y. 

Assumptions 

(i) i^ (i.e.; (si) G«). 

(ii) For every x G G^ - 1 , {y G G^ly ~ x} is finite. 

(iii) Points x, y G G*(°) for which there is a finite path x , xi, x 2 , . ■ ■ , x n 
with xo = x, x n = y, and (xjXj+i) G G^\ i = 0, . . . n — 1, are said to be 
connected 

(iv) We will assume that all connected components in G^ are infinite; 
or else that G^ is already connected. 

3 The Main Theorem 
3.1 The Graph Laplacian 

In this section we specify a fixed graph G (infinite in the non-trivial case) 
and an associated conductance function c. The associated graph Laplacian 
A c will typically be an unbounded Hermitian operator with dense domain. 

Our assumptions will be as above, and when the Hilbert spaces have been 
selected, our main theorem states that the graph Laplacian A c is essentially 
selfadjoint; i.e., the operator closure, also denoted A c , is a selfadjoint oper- 
ator. In sections 5-8 we obtain consequences and applications. 

The interpretation of this results in terms of boundary conditions will be 
given in section [7] below. It means that A c has a well defined and unique (up 
to equivalence) spectral resolution. Then the next objective is to find the 
spectrum of the operator A c . And a method for finding spectrum is based on 
"covariance." Covariance is used on other spectral problems in mathematical 
physics, and it offers useful ways of getting global formulas for spectrum. As 
we will see, infinite models typically have graph Laplacians with continuous 
spectrum. 

In the finite case, of course the spectrum is the set of roots in a charac- 
teristic polynomial, but unless there is some group action, it is difficult to 
solve for roots by "bare hands;" and even if we do, only the occurrence of 
groups offers insight. 
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A second approach to the finding spectra of graph Laplacians is "renor- 
malization:" Renormalization of hierarchical systems of electrical networks 
comes into play each time one passes to a new scale (upwards or down- 
wards). This requires additional structure, such as is found in iterated func- 
tion systems (IFSs), (see |BHS05j . |DJ07| . [JP98] . |Kig03| ), i.e., specified 
finite systems of affine transformations in Euclidean space that are then 
iterated recursively. 

When the mappings are so iterated on a given graph, the iterations 
may then be interpreted as scales in an infinite graph: (post-)composition 
of similarity mappings takes us further down the branches of a tree like 
structure in path space. We get martingale constructions as instances of 
renormalization. 

Theorem 3.1 The graph Laplacian A = is essentially selfadjoint. 

Proof. To get started we recall the setting. Given: 

G: a fixed infinite graph. (It may be finite, but in this case the 
conclusion follows from finite-dimensional linear algebra.) 



G= (G(°),G«). 



G^: the set of vertices in G. 



G^'i the set of edges in G. 



If x, y E G(°) is a given pair, we say that x ~ y when e = (xy) G G^ . 
For x £ G (0) , set 





Our standing assumptions are as follows: 

(a) nbh (x) is finite. 

(b) x nbh (x). 

U = l 2 (GW) = all functions v : G (0) -> C such that 




(3.2) 
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Set 

(u,v):= ^2 u(x)v(x), Vu,v e i 2 (g< >) . (3.3) 

:reG( } 

By H we refer to the completed Hilbert space £ 2 (G^). 

T>:= the set of all finitely supported v G H; i.e., t> is in ©iff 3F C G (0) , 
F = F v some /imte subset such that 

u(z) = 0, Vx G G (0) \F. 
e x '-= 8 X = Dirac mass, defined by 

c : G« - R + is a fixed function taking positive values. In network 
models, the function c is conductance; i.e., the reciprocal of resistance. 

Assumption (symmetry): c(xy) = c(yx) , V(xy) G G^. 

A = A (GjC) (3.5) 

is the Laplacian, and is defined on V as follows: 

(Av) (x): = Y,c (xy) (v (x) - v (y)) , G V, Vx G G^. (3.6) 

■ 

3.2 Lemmas 

We will need some lemmas: 

Lemma 3.2 The operator A is Hermitian symmetric on V, and it is posi- 
tive semidefinite. Specifically, the following two properties hold: 

{Au,v)p = (u,Av)(2, \/u,v G V; (3-7) 

and 

(u,Au} £ 2 > 0, VmG V. (3.8) 
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Proof. Both assertions are computations: 
In <ET~ 



(Au, v)(2 = J2 c ( x > y) \ u ( x ) - u (y)) v ( x ) 

= ^2 X) c ( xy ) u ( x ) v ( x ) ~ ^2 ^2 U ( y ^ c v ( x ) 

x^y ie~j/ 

= ^2 (x)c (xy) v (x) - ^2 (x)c (xy) v (y) 

x^y x^y 

= ^2 u ( x )^2 c ( x y)( v ( x ) - v (y)) 

xeG(°) y~ x 
= (u,Av) e 2. 

Note that the summation may be exchanged since, for each x € G^ ', 
the set of neighbors nbh (x) is finite. 
In dHSD, 

(u, Au) £2 

= ^2 5Z n ( x ) c ( u ^ ~ u ^ 

xr^y 

= J2 B c( x )\ u ( x )\ 2 -^2^2 u ( x ) c ( x y) u (y) > 

zeGM x~y 

where 

B c ( x ) = ^2c {xy) , x e G<® . (3.9) 

The second term in the computation may be estimated with the use of 
Cauchy-Schwarz as follows: Setting 

£ c = £ E c ^ \ u ^ - u ; ( 3 - 10 ) 

all x,y 
s.t. x^y 

we show that 

2{u,Au) =£ c {u) > 0. (3.11) 
Indeed using the conditions on c: — > M+ 

• c (xy) = c (yx) , V (xy) E G^ 1 ); 

• c(xx) = 0, Vx e GM; 
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• c(xy) > 0,V(xy) G G (1) , 
we get 

2(u,Au) = 2 £ c (x)Mx)| 2 -2^^^c(xy)u(y) 

= 2 ]T B c (x)\u(x)\ 2 -2Re^^^)c(xy)u(y) 



For the general theory of unbounded Hermitian operators and their ex- 
tensions, we refer the reader to |j0r78| , [Nel69j . |Sto51| . 

Definition 3.3 If A is an operator with dense domain D in a Hilbert space 
Tt, we define its adjoint operator A* by: 

A vector v is in the domain dom (A*) iff there is a constant K such that 

\(v, Au)\ < K \\u\\ , \/u G V. (3.12) 

When <\3. 12\ holds, then by Riesz, there is a unique w := A*v such that 

(w, u) = (v, Au), Vu G V. (3.13) 

Note that since T> is dense in Tt, w (=: A*v) is uniquely determined by 

Lemma 3.4 In the case of A = A/qa an d Tt = I 2 (G^), the vector A*v 
for v G dom (A*) is given by the expression 

(A*v) (x) = Y c (?y) (« (*) - v (y)) . (3.14) 

Proof. Since the sum in (|3.13p is finite, the RHS is well defined if v G 
dom (A*). Since A*v G Tt, 

Y |(A^)(x)| 2 <oo. (3.15) 
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Set w (x):= Yly yx c [ x v) ( v ( x ) - v (y))- 

We claim that (|3,12p then holds. Indeed 



Yl {^2 c ( x y)( v ( x )- v (y))) u ( x ) 

;e<3(°) \y~x ) 



{w,u)p = 

= ^2 v ( x )J2 c i x v) ( u ( x ) - u (y)) 

x y^x 

(by the exchange of summation and Lemma 13. 2p . ■ 

Lemma 3.5 Let (G, c) and A = A(£ jC ) here as in the previous lemma. 
Then the equation 

Av = -v (3.16) 
does not have non-zero solutions v £ £ 2 (G^) . 

Proof. It is immediate from (13. 71) in Lemma 13.21 that eq. (|3.16j) does not 
have non-zero solutions in T>, but the assertion is that there are no non-zero 
solutions in any bigger subspace. 

Also note that every solution in £ 2 (G^ )) to eq. (|3.16p must be in 
dom (A*), i.e., the domain of the adjoint of A with D as domain. 

If v : -> C is a solution to (1536]) . then 

v (x)Av (x) = — \v (x)\ 2 , Vx £ G (0) (3.17) 

which yields v(x)Av (x) < 0, Vx G Hence £ c (v) < 0; see (^T0l) - (^TTj) . 

But by (|3TT1) . then £ c (v) = 0. 

It follows from (|3.10p that v must be constant on every connected com- 
ponent in Since all the connected components are infinite, v must be 
zero. ■ 

Remark 3.6 We stress that (jg.iffi) may have non-zero solutions not in £ 2 . 
For these solutions v, the energy will be unbounded. 

Example 3.7 Let a graph system (G,c) be determined as follows: 

G(°) = No = {0,1,2,...}, 

G« : nbh(0) = {1}, 
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nbh (n) = {n — 1, n + 1} if n > 0, and 
c (n, n + 1) = n + 1. 
T/ien i/ie Laplace operator A c ifiZZ 6e unbounded in £ 2 as follows from 

\ 



1 


-1 











-1 


3 


-2 











-2 


5 


-3 











-3 


7 


-4 



-n 2ra + 1 - (n + 1) 



T/ten 



(An) = wo — «i, and 

(Au) n = (2n + 1) u n - nu n -i - (n + 1) Vra > 1. 

For solving (|3.16p . initialize vq = 1. Then 

vi = 2^o = 2, 

W2 = -, and inductively 



u n+ i = 2v n - 
We get vi < v 2 < ■ ■ ■ < f n -i < v n < 



v n +l > 2 



n + 1 
• and 

n 

n + 1 



v n -i- 



Hence for the truncated summations for I 2 and £ applied to this solution v ; 
we get 



1 * 

-S AT (V) = - ^ < 



-JV 



fc=0 



which tends to — oo. 

The following lemma is from the general theory of unbounded operators 
in Hilbert space, [Nel69j . |Sto51j . [vN3T] . 
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Lemma 3.8 Let A be a linear operator in a Hilbert space 7i and defined in 
a dense domain T>. 

Then A is essentially selfadjoint (i.e., has selfadjoint closure A) if the 
following conditions hold: 

(i) (u,Au) > 0,Vu G P 

(ii) dim{u G dom (A*) \A*v = —v} = 0. 

Proof. This is in the literature, e.g. [vN31| . The idea is the following, if (i) 
is assumed, then there is a well defined bounded operator 

T=(l + Ay 1 

precisely when (ii) is satisfied. ■ 

In our analysis of the graph Laplacian A c in fl3.6f) we shall need one 
more: 

Lemma 3.9 Let A c be as in Then for all v £ T>, 

(A c v)(x) = 0. (3.18) 

Ln fact, when v is fixed, the number of non-zero terms in (\3.18\i is finite. 

Proof. The finiteness claim follows from the assumptions on (G, c) we listed 
in section [5J 

A direct computation yields the result: 

(Acu)(a;) =^2^2c(xy)(v(x)-v(y)) 

x y xr^y 

= J2v{x) B c (x)-^2v (y) c (v x ) 

x y x^y 

= Y,v{x)B c (x)-Y,v{y)B c (y) 

x y 

= 0. 
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4 Operator Theory 



Once the operator theoretic tools are introduced, we show in section [5] below 
that class of infinite graph systems (G, c) where G is a graph and c is a 
conductance function (the pair (G, c) satisfying the usual axioms as before), 
have the spectrum of the associated Laplace operator A c continuous. This 
refers to the £ space of G(°), i.e., the Hilbert space is l 2 (G^) where as 
usual G^ denotes the set of vertices. 

It is important that G^ is infinite. Otherwise of course the spectrum is 
just the finite set of zeros of the characteristic polynomial. See Example 16.61 
below. 

We give an operator theory/spectral theory analysis, with applications, 
of a class of graph Laplacians; and we have been motivated by a pioneering 
paper jPow76] which in an exciting way applies graphs and resistor networks 
to a problem in quantum statistical mechanics. In one of our results we 
establish the essential selfadjointness of a large class of graph Laplacians on 
graphs of infinite networks. (A Hermitian symmetric operator with dense 
domain in Hilbert space is said to be essentially selfadjoint if its closure is 
selfadjoint, i.e., if the deficiency indices are (0,0). See Definition 14. II below! 
There are many benefits from having the graph Laplacian A essentially 
selfadjoint. 

Here is a partial list: 

(a) We get the benefit of having the spectral resolution for the selfadjoint 
closure, also denoted A for notational simplicity. 

(b) We get a spectral representation realization of the operator A, i.e., 
a unitarily equivalent form of A in which an equivalent operator A~ may 
occur in applications. See e.g., [Arv02j . [PS72J. 

(c) We get a scale of Hilbert spaces, TL S for s in R, defined from the graph 
of the operator A s where the fractional power A s is defined by functional 
calculus applied to the selfadjoint realization of A. See [Jor04j. 

(d) Gives us a way of computing scales of resistance metrics on electrical 
networks realized on infinite graphs, extending tools available previously 
only for finite graphs; see [BD49J. 

(e) The case s = 1/2 yields an exact representation of the energy Hilbert 
space associated with a particular system (G, c) and the corresponding graph 
Laplacian A = A( G c ). 

(f) Gives us a way of computing fractional Brownian motion on graphs, 
allowing an analytic continuation in the parameter s, and with s = 1/2 
corresponding to the standard Brownian motion; see e.g., [DJ07J . [Jor06j . 

In the course of the proofs of our main results, we are making use of 
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tools from the theory of unbounded operators in Hilbert space: von Neu- 
mann's deficiency indices, operator closure, graphs of operators, operator 
domains, operator adjoints; and extensions of Hermitian operators with a 
dense domain in a fixed complex Hilbert space. Our favorite references 
for this material include: [AC04j . U0r77| , |j0r78| , [JPOO] . |Nel69j . [vN3T] . 
[Sto51| . For analysis on infinite graphs and on fractals, see e.g., [BHS05J, 
[CS07] . [DJ06] . |HKK02j . |Hut81j . [JP98] . [JKS07| . |Kig03| , |BD49j . 

Definition 4.1 Let A be a Hermitian linear operator with dense domain T> 
in a complex Hilbert space 7i. Set 



where i = y— T. Then the two numbers n± := dimP-t are called the defi- 

ciency indices. 

Von-Neumann's theorem states that the initial operator A is essentially 
selfadjoint on T> if and only if n + = ra_ = 0. It has selfadjoint extensions 
defined on a larger domain in H if and only if n + = n_. 

The following two conditions on a Hermitian operator, (A) and (B), 
individually imply equal deficiency indices, i.e., n + = 

(A) For all v £ T>, we have the estimate 



i.e., A is semibounded. 

(B) There is an operator J : 7i — > 7i satisfying the following four condi- 
tions: 

(i) J (u+ av) = Ju + aJv, for Vu, v £ 7i, a £ C 

(ii) (Ju, Jv) = (v,u) \/u,v £ Tt; (J is called a conjugation!) 

(iii) J maps the subspace V into itself, and 



(iv) J 2 = id; J is of period 2. 

Remark 4.2 There are many examples (see the Appendix) where either 
(A) or (B) is satisfied but where the operator A is not essentially selfadjoint. 



V± : = {v± £ dom (A*) | A*v± = ±iv±} 




(v,Av) > 0, 



JAv = AJv, \fv £ V. 
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Both conditions (A) and (B) hold for a graph Laplacians A c , and The- 
orem [3J] states that A c is essentially selfadjoint. 

For Riemannian manifolds with boundary, there is a close analogue of 
the graph Laplacian A c above; but it is known (see section [7J) that these 
continuous variants are typically not essentially selfadjoint. 

Indeed the obstruction to essential selfadjointness in these cases captures 
a physical essence of the metric geometry behind the use of Laplace operators 
in Riemannian geometry. 

5 The Energy Form 
5.1 Operators 

In section [3] we proved essential selfadjointness of the graph Laplacians A c . 
This refers to sequence space £ 2 , the Hilbert space of all square-summable 
sequences indexed by the points in G^°\ and equipped with the usual £ 2 - 
inner product. 

This means that the axioms for A c are such that boundary conditions at 
infinity in G are determined by computations on finite subsets of the vertices 
in G. (In the Appendix, we will contrast this state of affairs with related 
but different boundary conditions from quantum mechanics.) Recall that 
A c is generally a densely defined Hermitian and unbounded operator. So 
in principle there might be non-trivial obstructions to selfadjointness (other 
than simply taking operator closure.) Recall (Definition 14. ip that a given 
Hermitian operator with dense domain is essentially selfadjoint if and only 
if the dimension of the each of the two "defect eigenspaces" is zero. 

So that is why we look at the "minus 1 eigenspaces" for the adjoint 
operator, A*u = —u. 

For potential theoretic computations we need an additional Hilbert space, 
the Energy Hilbert space Hg (details below.) For example the voltage po- 
tentials associated with a fixed graph Laplacian are typically not in £ 2 (G^) 
but rather in an associated Energy Hilbert space. Our Laplace operator 
A is formally Hermitian in both the Hilbert spaces £ 2 and H.£(the energy 
Hilbert space). We show that the Laplace operator A is essentially selfad- 
joint both in I 2 and in Tig. In both cases, we take for dense domain T> the 
linear subspace of all finitely supported functions G^ — > R. 

Our setting and results in this section are motivated by [Pow76] and 
|BD49j . 

There are several distinctions between the two Hilbert spaces: For ex- 
ample, the Dirac functions {5 x \x G G^} form an orthonormal basis (ONB) 
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in £ 2 , but not in Tig. The implication of this is that our graph Laplacians 
have different matrix representations in the two Hilbert spaces. In speaking 
of "matrix representation" for an operator in a Hilbert space, we will always 
be referring to a chosen ONB. 

We shall need the operator A in both guises. One reason for this is that 
for infinite graphs, typically the potential function v solving Av = 5 X — 8 y , 
for pairs of vertices will not be in I 2 , but nonetheless v will have finite 
energy, i.e., S (y) < oo, meaning that the energy form applied to v is finite. 
Caution: The sequence v might not be in the ^ 2 -space. Specifics in Example 
ICT below! 

When we study the Laplace operator A, our questions concern its spec- 
trum, and its spectral resolution. The spectrum will be contained in the 
half- line [0, oo), but (as we show in examples) it can be unbounded, and it 
can have continuous parts mixed in with discrete parts. In case the conduc- 
tance function is "very unbounded," as an operator in l 2 , it may be necessary 
to pass to a proper operator extension in an enlarged Hilbert space to get a 
different selfadjoint realization of A. 

These operator issues only enter in case A is unbounded. The unbound- 
edness of the operator A is tied in closely with unboundedness of the con- 
ductance function c which is used. Recall A = A c depends on the choice of 
c. If c is unbounded "at infinity" (on the set G^ of edges), then the resis- 
tors tend to zero at distant edges. Intuitively, this means that "the current 
escapes to infinity," and we make this precise in the language of operator 
theory. Our general setup here will be as in [Pow76]. 

We will need a second Hilbert space, the energy Hilbert space Tig. Here 
the inner product is the energy quadratic form. Since the energy form 
evaluated on a function v is defined in terms of the square of differences 
v(x) — v(y), it follows that the elements in Tig are really sequences modulo 
constants. 

Let G = (G^°\ G^') be a graph satisfying the axioms in section [3l with 
vertices and edges . Let 

c : G (1) -> M+ (5.1) 

be a fixed function (called conductance.) If e = (xy) £ G^ 1 ' , we say that x ~ 
y, and the function c must satisfy c(xy) = c(yx), symmetry. In particular, 
for a pair of vertices x,y, c(xy) is only defined if x ~ y, i.e., if (xy) £ G^\ 
For every x £ G^°\ we assume that 

nbh(x) = {y £ G (0) |y ~ x} (5.2) 
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is finite, and that x nbh(x). 

Following eq. (|3.9p . we study functions v : — > C for which 

£c{v)= ^2 c(xy)\v (x) - v (y)\ 2 < oo. (5.3) 

all x,y 
s.t. x^y 

Clearly we must work with functions on G^ modulo constants. Setting 

£ c (u, v) ■. = ^2 c ( x v) ( u ( x ) - u (y)) ( v ( x ) - v (y)) , ( 5 - 4 ) 



all x,y 
x^y 



we get 



\£ c (u,v)\ 2 < £ c {u)E c {v) , Mu,veV, (5.5) 



by Schwarz' inequality. 
Setting 

(u,v)e : = £ c (u,v) (5.6) 

we get an inner-product, and an associated Hilbert space Tig of all functions 
v for which (|5.3p holds. 
The triangle inequality 

S c (u + v)5 <S c (u)5 +g c (v)5 (5.7) 

holds; or equivalently 

\\u + v\\ £ < \\u\\g + \\v\\ £ , V-u, v G Tie- 

In the next result we give a characterization of the Hilbert space 7i.£ directly 
in terms of the selfadjoint operator A c from sectional 

Recall A c is the closure of the operator A c with dense domain T> in 
i 2 (G^Y It will be convenient to write simply A c for the closure. Since it is 
selfadjoint, we have a Borel functional calculus; i.e., if / is a Borel function 
on R, and if P (•) is a projection valued measure for A c , then 



A c = / XP(dX), (5i 

and 



/(A c ):= J f(X)P(d\). (5.9) 
For the corresponding (dense) domains, we have 

|A| 2 \\P{dX)vf < ool , (5.10) 



dom(A c ) = jt> G t 2 ( Gi0) )\£ 



17 



and 



dom (/ (Ac)) = L e i 2 (GW) I 1/ (A)| 2 ||P (dA) vf < ooj . (5.11) 

Theorem 5.1 Let G = ((jM^W) and c : GW -» M+ 6e as described, and 
let Tig be the energy Hilbert space. Let A c be the selfadjoint graph Laplacian 
in t 2 (G<°)) /roro Section H 
(a) T/ien 

dom(Ay 2 )=^n^ 2 (G(°)). (5.12) 



(6) In general the right hand side is a proper subspace of I 2 (G^ )). 
Proof. We proved in (|3.1ip . Lemma [321 that 

S e (u,v) = 2(u,A c v) (5.13) 
for all functions u, v on G^ for which the two sides in (|5.13p are finite. 



A)J\ 2 



iff v € 



If u = v, then the expression on the RHS in f)5. 13j) is 2 
dom ^A^ 2 ^ . This follows from (|5.1ip applied to 

/ (A) : = V\, A G [0, oo). 

Since dom(A c ) C dom (aJ /2 ) C I 2 (G^) the desired conclusion (pU2|) 
holds. 

To see that dom [A l J 2 ] may be a proper subspace of 



H £ = {v\£ c (v)<oo}, (5.14) 

consider the following example (G, c) built on the simplest infinite graph 
G(°) : = Z. ■ 

Example 5.2 

G(°) = Z, 

GW = {(n,n±l)|neZ} ) ana 
c:GW^R + , c=l onG^. 

The corresponding graph Laplacian is 

(Av) (n) = 2v (n) - u (n - 1) - v (n + 1) , Vn € Z. (5.15) 

If k & Z + is given, we claim that there is a unique function v £ solving 

Av = S -8 k . (5.16) 
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Existence: Set 

{0 if n < 
-n ifO <n<k (5.17) 
—k ii k <n. 

A substitution shows that the function v in (|5.17|) satisfies (|5.16|) . 
Uniqueness: Let w £ TLs be a solution to (|5.16p . Then 

£ c (u, v — w) = 0, VuGP. 

Since G = (G^°\ G'- 1 )) is connected, V is dense in Kg; and so the difference 
v — w must be a constant function. But the Hilbert space TLs is defined by 
moding out with the constants. Hence, v = w in TLs- 
The following three facts follow directly from (j5.17|) : 

(i) v is non-constant; 

(ii) v £l 2 {1); and 
(hi) £ c (v) < oo. 

In fact, an application of (|5.13p yields 

£ c (v) = 2k. (5.18) 

Proof of flSTTB]) : 

£ c («) = 2(u,Aw) by (I5TT31) 
= 2(u,5 - £*) 
= 2(v(0) — v (Jb)) 
= 2A; by (j5TT|> . 

Theorem 5.3 Zei G = (G (0) ,G (1) ) and c : G (1) -> R + be a graph system 
as in the previous theorem, and in section^ i.e., we assume that the pair 
(G, c) satisfies the axioms listed there. Let A c be the corresponding graph 
Laplacian with a choice c for conductance. 

Let a, j3 € G^ be a fixed pair of vertices. Then there is a unique function 
v £ TLs, i.e., £ c (v) < oo satisfying 

A c v = 5 a -5 p . (5.19) 

The solution to is called a voltage potential. Moreover, 

£ c (v)=2{v{a)-v((3)). (5.20) 
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Proof. The argument for uniqueness is the same as in the previous proof. 

To prove existence, we will appeal to Riesz' theorem for the energy 
Hilbert space Tig. 

Hence, we must show that there is a finite constant K such that 

\u (a) - u (f3)\ < KE C («)* for all u G V. (5.21) 

Motivated by Ohm's law, we set Cl(e) := c(e) _1 , Ve G £?W. By the 
assumptions in section [3J we may pick a finite subset xo,xi,X2, . . . ,x n in 
G<® such 

x = a, x n = (3, and , _ 

e { = (xi x i+ i) G G (1) , i = 0, 1, . . . , n - 1. 

Then 

n-l 

a) - it | < ^ |u (x^ - u (xj+i)| 



\ u , - 

i=0 



i i 

(n— 1 \ 2 /n— 1 \ 2 

^ ^ (ej) J {^2c(ei)\u(xi)-u(x i+ i)\ 2 \ (by Schwarz) 

(n-l \ 2 

To get a finite constant -fT in (|5.2ip . we may take the infimum over all 
paths subject to conditions (|5.22j) . connecting a to (5. 

An application of Riesz' lemma to Tig yields a unique v G Tig such that 
for all u G T>, we have the following identity: 

it (a) -u(/3) = -S c (y,u) 

= (A c v,u) (by (E331)). 
Using again density of T> in Tig, we get the desired conclusion 

A c v = 5 a -5(3. (5.23) 



Corollary 5.4 Let (G,c) satisfy the conditions in the theorem. Let a,f3 G 
G( ) ; and Zei u G 6e i/te solution (potential) to 

A c v = S a - dp. 
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Then 

n-l 

£ c (v) < 2 inf V H (e») (5.24) 

where eo, ei, . . . , e n _i G G^ 1 ) is a system of edges connecting a to (3, i.e., 
satisfying the conditions listed in 

Proof. This follows from the previous proof combined with the fact that 

sup \£ c (u, v)\ 2 = £ c (v) . (5.25) 

£c{u) = l 



5.2 A matrix representation 

While A c may be understood as an operator, it is also an oo x oo matrix. 
Since the set nbh (x) C is finite for all x G G^°> , A c is a banded matrix. 
To see this, note that when x G G^ is fixed, the summation 

(A c v) (x) = Y,c (xy) (v (x) - v (y)) (5.26) 

y~x 

is finite for all functions v : G^ — * C. 

Since G is assumed connected, the only bounded solutions v to the equa- 
tion 

A c v = (5.27) 

are the constants. 

Solutions v to (I5.27P are called harmonic, or c-harmonic. 

There are examples of systems (G, c) which are connected and have un- 
bounded non-constant harmonic functions, e.g., models with G^ = Z 3 , or 
tree-models. 

In the general case, introducing 

B c {x) :=Y,c(xy), x£ G®; (5.28) 

we see that (|5.26p takes the following form 

(A c v) (x) = B c (x) v(x)-J2 c fry) v (y) ■ (5-29) 
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Hence eq. (|5.27p may be rewritten as 



v ( x ) = -^-r^^2 c ( x y) v (y) ■ ( 5 - 3 °) 

0C W y~ X 

It follows that harmonic functions on G^ satisfy a mean value property. At 
every x € Gf(°) formula (O0|) expresses i> (x) as a convex combination of its 
values on the set nbh (a;). 

In matrix language, x — > Z3 C (x) represents the diagonal matrix-entries; 
and c(xy) the off-diagonal entries. Since {y € G^\c(xy) ^ 0} is finite, we 
say that the matrix for A c is banded. It is clear that products of banded 
matrices are again banded; and in particular that the summations involved 
in matrix-products of banded matrices are all finite. Hence, each of the 
operators A c , A 2 , A^, . . ., is banded. Since by Theorem 13 .l\ A c is selfadjoint 
as an operator in I 2 (G (0) ), the fractional power A^ 2 is well defined by the 
Spectral Theorem. The matrix-entries of A^ 2 are the numbers 

(x) , x,y € G^. (5.31) 

1/2 

It can be checked that if G is infinite, the matrix for A c ' is typically not 
banded. The same conclusion applies to A* when s £ M\N. 

5.3 Example 5.2 revisited 

The system (G, c) in Example 15.21 does not have non-constant harmonic 
functions. This can be seen from the representation of A (in Ex. 15. 2p as a 
ZxZ double infinite matrix, i.e., 

(Av) (n) = 2v (n) — v (n — 1) — v (n + 1) 

= v (n) — v (n — 1) + v (n) — v (n + 1) 

= ^2 v ( n ) — v ( m ) , n £ Z. 
In matrix form, A from Example 15.21 is as follows: 



22 



-1 












-1 

2 
-1 










-1 

2 
-1 









-1 

2 
-1 



V 



Using Fourier series 
f{x) 



J2 v ( n ) einx e ^(-^o 



I / 0*01 dx; 



we arrive at the representation 

'A/) (a?) = 2(1 -cosx)/(x) 
= 4sin 2 (|)/(x), 

proving that A has Lebesgue spectrum, and 



spec (A) = spec ( A 

I 2 L 2 



5.4 Banded Matrices (A Preview) 



[0,4]. 



(5.32) 
(5.33) 

(5.34) 



(5.35) 



It is immediate from the matrix representation for A c in Example 15.21 that 
it has a banded form. We will take up banded infinite matrices in detail in 
section [8] below. 

1/2 

Since A c is selfadjoint, its square-root A c is a well defined operator. 
However its matrix representation is typically not banded; see (I5.3ip . For 
A c ' in Ex. 15.21 one can check that the (m, n)-matrix entries are 



4 • (n — m) + 1 
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5.5 Extended Hilbert Spaces 

To understand solutions v to operator equations like 

A c v = 5 a - 5(3 

as in (|5.23p . potential functions it is convenient to extend the Hilbert space 
£ 2 (GW). Indeed we saw in Example 15.21 that the solutions v to equations 
like (103]) are typically not in £ 2 (G<W). 

Definition 5.5 The space TC C (s). 

A function v : — > C is said to belong to the space 7i c (s) if there is 
a finite constant K = K (s) such that the following estimate holds: 



£ v (x) (A s c u) (x) 

;eG(°) 



< -RT(s) ^ |u (x)| 2 for all «eP. (5.36) 



If ()5.36p holds, then by Riesz, there is a unique w £ £ 2 (G^) such that 
v(xj (A'n) (x) = for all « G V ( C ^ 2 ); (5.37) 

and we set 

IMI W M : = Ikll^(GO) • (5.38) 

By abuse of notation, we will write A^u = w when v S TC C (s). 

If two functions vi for z = 1, 2 are in W c (s), and if A^Vi = Wi £ t 2 (G^) , 
we set 

("1^2)^(8) : = (wi,w 2 ) e 2 (5.39) 
= ^ wi {x)w 2 (x) . 

Remark 5.6 We proved in section\3\that 

^2 (A c n) (x) = for Vu G P. 

Hence the constant function v\ (x) = 1 on G^ is in 7i c (1), and H^iH^m = 
0. Hence in considering the extension spaces, we shall work modulo the 
constant functions on G^ . 
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Theorem 5.7 For every s £ K, the space 7i c (s) is a Hilbert space when 
equipped with the inner product (|5. 39\i , and the norm (|5, 38\i . 

Proof. The idea in the proof follows closely the construction of Sobolev 
spaces, by analogy to the continuous case. The key step in the verification 
of completeness of 7i c (s) is the essential selfadjointness of A c as an operator 
m I 2 As before, we use the same notation A c for the closure of A c , 

defined initially only on the subspace T> in £ 2 (G^). Formulas (|5.9p - (|5.1ip 
above now allow us to define the selfadjoint operator A*; and this operator 
is closed in the sense that its graph is closed in I 2 x I 2 (G(°)). The 

completeness of 7i c (s) now follows from this, and an application of Riesz; 
see the estimate (|5.37|) . ■ 

Corollary 5.8 Let (G, c) be an infinite graph, and let c : — > be a 

conductance function satisfying the axioms above. Let a,/3 S G^°\ and let 
v : G^ — > C be a solution to A c v = S a — 6/3; i.e., to <{5.23 . 
Assume v £ Tig. Then 



veHc 



and we have 



and 



1 1/2 



nWc(i); 



(5.40) 



(5.41) 



(5.42) 



Proof. To prove (|5.40p . we must check the a priori estimate (|5.36p for 
s = 1/2, and s = 1: 



Verification of rfBTM]) for s = 1/2 

Let v satisfy the stated conditions, and let u £ T>. Then 



v(x)(Aiu)(x) 



v(x)A c A c 2 U (x) 

;eG(°) 



£ c ( v, A c 2 u 



1 1 / -1 \3 

< -£ c (v) 2 £ c I A c 2 u j (Schwarz) 

1 c f \± 11 11 

= -^ £ W 2 IMIf2( G (0)) , 
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where we used the identity 



£ r [ A r 2 u) =2 lilt 



l<2(GC°)) 

2 £ |u(x)| 2 

rcGG(°) 



valid for \/u £ V. 



Verification of (I5.36|) for s = 1 

With v and u as before, we must estimate the summation: 



^ v (x) (A c it) (x) 

c eG(°> 



y~] (A c u) (x) u (x) 



^2 ($a 0) - 6p (x)) u (a 

XGG(°) 

= |-u (a) — -u (/3)| 

< 2 ||n||^( G (o)) , VueD. 



Once (|5.40p has been checked, the exact formulas (|5.4ip and (15.421) follow: 
Firstly, 



and secondly 



\v l 

2 



Aft, 



£ 2 

1 1 



|^|| ! = \\A c v\ 
= \\5 n — 5 



a — 
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Remark 5.9 In conclusion ( |5.^0fl in Corollary \5.8\ is not best possible. In 
fact, the optimal range of the fraction s for which the potentials v are in 
H. c (s) may be computed explicitly in Example \5.S\ and related examples. De- 
tails in the next subsection. 

In Example\SM G (0) = Z, = {(n,n± 1) \n £ Z}, and c = 1. Let 
k £ N. The graph Laplacian A is given in formula (|5, 15\i . 

Let v be the unique solution to the potential equation 

Av = 5 -5 k . (5.43) 
Then v £ Tt(s) if and only if s > 1/4. 
Proof. Setting 

v(z) = v nZ n , an d z = e lx , (5.44) 



we get 



and therefore 



z (z k — l) 

v(z)= f J ; (5.45) 
(z- 1) 



Ism ( -g- ] 



\v(x)\ = \ 2 )j/ n [ . (5.46) 



k 2 > 

sva 1 (x/2) 

Since, in the spectral representation, the graph Laplacian A is multiplication 
by 4 sin 2 (x/2), the question: "For what exponents s is 

v£H{s)T (5.47) 

is decided by the asymptotics near x = of the function (A s v) (x). Using 
(|5.46p . we see that A s v is in I? (— vr, n) if and only if x 2s_1 £ L 2 near x = 0; 
and this hold if and only if 

s > -7 (5-48) 
4 

as claimed. ■ 



5.6 Lattice Models 

Example 5.10 We proved that potential functions are often not in I 2 (G^ ^) , 
but in general the problem is more subtle. 
The setting is a follows: 

G= (G(°),G«) a given graph; 
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c : G(i) _> R+ 

a ^iwen conductance function; 
A c = the corresponding graph Laplacian; 
a, (3 € G^ ) a /ixed pair 0/ vertices, a ^ (5. 

With this, we say that a function v : G^ ^ — > R zs a potential if 

A c v = S a -S p . (5.49) 

In the next result we show that lattice models iP with D > 2 have £ 2 
potentials. 

5.7 Preliminaries 

By 2/° we mean the rank D-lattice of vertex points n = (ni,n,2, ■ ■ ■ ,tid), 
ni 6 Z, i = 1,2, ...,D. Every point n G lP has 2D distinct nearest 
neighbors 

(ni, . . . , ni ± 1, n i+ i, . . . , no) , (5.50) 

so nbh (n) consists of these 2D points; and is the corresponding set of 
edges. In the discussion below, we pick the constant conductance c = 1, i.e., 
a system of unit-resistors arranged in nearest-neighbor configurations. See 
Fig. 1 for an illustration of the simplest lattice configuration, D = 1,2, and 
3. 



Fig. la: D = 1 
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Fig. lb: D = 2 



Fig. lc: D = 3 

Fig. 1. Lattice configurations in the rank-D lattices 7lP 
with nearest-neighbor resistors. 

Proposition 5.11 The potential functions v, i.e., solutions to fl5,^ff[ ) with 
c = 1 are in £ 2 (Z D ) if D > 2. 

Proof. Recall that the .D-torus T D is the compact dual of the rank-D lattice. 
Pick coordinates in T s.t. x = (x\, . . . ,xd), — tt < Xi < tt, % = 1, 2, . . . , D. 
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Then, by Parseval, 

e (Z D ) ~ L 2 (T D ) . 

By the argument from Example [521 we see that A has the following spectral 
representation in L 2 ((— 7r,7r] x • • • x (— 7r, 7r]) 



D times 

D 



(Av)(x) =4^sin 2 jv(x). (5.51) 

k=l 

Introducing spherical coordinates in R D , we get the representation 

dx = p D - 1 dS l (5.52) 

where 

and whence dSi denotes the rotationally invariant measure on the sphere in 
R D . 

The question of deciding when the solution v to (I5.49|) is in t 2 {IP) can 
be better understood in the spectral representation v (x) for x = (xi, ■ ■ ■ , xd) 
close to 0, i.e., p{x) ~ 0. 

Using (|5.5ip - (|5.52p . we see that the potential function v is in £ 2 if D > 2. 

More generally, the argument from Example 15.21 proves that in iP , the 
potential function v is in H (s) if s > ^p-- ■ 

The next results illustrate new issues entering the analysis of Z^-graphs 
when D > 1, compared to the D = 1 case. 

Corollary 5.12 For the case D = 3 in the lattice model in Example 15. 10[ 
consider k = (ki, &2, ^3) £ Z 3 \ (0) fixed, and let 

v k : 1? -» M 

the solution to the potential equation 

Av k = 5o-5 k . (5.53) 

Then 

lim v k (n) = 0. (5.54) 
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Proof. Our notation is as follows: n = (ni,n2,ns) G Z 3 and by "n — > oo" 
we mean: 

\n\ = ^n\ + n\ + n| — > oo. 

Moreover, 

5 k : I? -» R 

is the usual Dirac mass 

<5 fc ( n ) = ^fc^m^j.ni^fcs.ns- (5-55) 
We proved in Proposition 15.111 (D = 3) that 

^ K (n)| 2 f= IKI| 2 ) < oo; (5.56) 

and so in particular, the conclusion (|5.54p must hold. ■ 

Our next example illustrates that the potential equation (|5.53|) has un- 
bounded solutions in case D > 3. This will also provide concrete cases 
of unbounded harmonic functions, i.e., functions w : lP — ► R for which 
Aw = 0. 

To aid the construction, we include the following lemma which is about 
the general case of systems (G, c) as analyzed in sections [3][H above. 

Lemma 5.13 Let G = (G {0 \G^), and c : GW -> R+, be a graph system 
as described in Theorem \5.1\ and let A c be the graph Laplacian. 

Let a, (3 G G^ be given, a ^ [3. Then there is a 1-1 correspondence 
between two classes of functions v : G^ — > R, and functions I : G^ — * R, 
where the two classes are given as follows: 

Class 1. 

A c v = 5 a -5f3 (5.57) 
Class 2. I : G^ — > R satisfying (Kirchoff's Laws): 

(«) Z y ~ x I (x, y) = (5 a - 6/3) (x) , Vx G G<°) ; and 

( & ) Si c (^i^+i)" 1 -f (aw+i) = /or a// x Q ,x 1 ,X2, ... ,x n £ G (0) stifr- 
jeci to xo = x n , and X{ ~ Xi+i, i.e., all closed loops in G^ . 

The connection between the two classes is given by the following formula: 

c (xy) -1 / (xy) = v (x) — v (y) , V (xy) G G^ . (5.58) 
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The function v is determined from I uniquely, up to a constant, when I 
is known to satisfy (a) -(b). Moreover, 



Proof. Left to the reader. The arguments are included in the proof of 
Theorem 15.11 ■ 

Example 5.14 The T> = 2 lattice model; i.e., G^ = I?, edges given by 
nearest neighbors as in Fig. lb; and c = 1. 

We consider the equation (|5.57p for a = (0,0) and (3 = (1, 1). The two 
different solutions v to (|5.57p will be presented in the form of Class 2 in 
Lemma 15.131 i.e., in terms of current functions defined on the edges in G. 
First recall that the Laplace operator A in the Z 2 -model is 



We now describe the two current functions / which correspond to the 
two solutions to (I5.60|) . 



e 




(5.59) 



(Av) (m, n) = 4v (m, n) — v (m — 1, n) — v (m + 1, n) 
—v (m, n — 1) — v (m, n + 1) , V (m, n) £ 1? . 



Eq. (|5.56p then takes the form 



(5.60) 
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> > > > < < < 

A A A A V V V 

> > > > < < < 
AAA A VVVV 

> > > > < < < 
AAA A VVVV 

< < < > < < < 
AAA A VVVV 

< < < > (i,i) < < < 

A A A A A A A A 

< < < (0,0) > > > > 

VVVV V AAA 

< < < > > > > 
VVVV V AAA 

< < < > > > > 
VVVV V AAA 

< < < > > > > 



Fig. 2. The function I for the first solution v to (|5.60p . 

Flow design for the current function /. The symbols ">" indicate 
arrows in the direction of the current flow. An arrow points in the 
direction of voltage drop. 
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And now the (different) function I for the second solution to f|5.60|) : 



t o \ 

v 



< 2 < 2 

2 2 2 

v v 

<\ <\ 

\ o i i 

v v 

< \ <\ <\ 

A 

\ i i 

4 4 2 u u 

V V 

1 > (i,i) < \ 

A 

1 1 1 

2 4 4 U U 

A V 

U o,o) > ^ j > ^> ^> ^> tj> 







2 4 ^ 2 ^ 2 ^ 2 ^ 2 

••• 



Fig. 3. 

Flow design for the current function /. The symbols ">" indicate 
arrows in the direction of the current flow. An arrow points in the 
direction of voltage drop. 

5.8 The Resistance Metric 

Let G = (GW,GW) be a graph satisfying the axioms from section O and 
let 

c : -> R+ 

be a conductance function. Let £ c (•) be the corresponding energy form, and 
let A c be the graph Laplacian. 

Pick a definite point in the vertex set G^ '. Now for every x £ G^ let 
v x S TLs c be the solution to 

A c v x = 5o-5 x . (5.61) 
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Set 

dist (x,y) : = £ c (v x - v y )i (5.62) 



| Vx Vy 1 1 £ 



for x, y G G^ ^. We say that x,y — > dist c (x,y) is the resistance metric on 
G*- ). It is immediate from ()5.61f) that it satisfies the triangle inequality. 

Proposition 5.15 T/ie following formula holds for the resistance metric: 

di st( X ,y) = V2(v x (y) + Vy ( X )-v x ( X )-v y (y^. 

Proof. In view of (|5,62p . it is enough to compute E c (v x — v y ) for pairs of 
points x, y. 

Let x, y e be given, and let v x ,v y be the potential functions from 
flESH) . Then 

£ C («a: - = 2 (A C ("Vx - V y ) , V-,. - Uj,)^ 

= 2(5o -S x — (8q - 5y) , v x - v y )p 

= 2 {{v x - v y ) (y) - (v x - v y ) (x)) 

= 2 (v x (y) + v y (x) - v x (x) - v y (y)) . 



Example 5.16 (See also Example \3. 7[ ) Let A c be given by the following 
00 x 00 matrix: 



( 



1-1 

15 -2 2 
_2 2 2 2 + 3 2 -3 2 







-n 2 n 2 + (n+l) 2 -(n+lf 







So G(°) = N , G (1) = {(0, 1) , • • • , (n - l,n) , (n, n + 1) , • • • }, and c (n, n + 1) 
(n + l) 2 . T/ie /irsi vertex has one neighbor, and the later two. 
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The potential equation ( 15.62j) may be solved by inspection, and we get 
the following formula for the resistance metre dist c in Proposition I5.15t If 
m < n (in No) then 

,. / N ( 1 1 

dist c (m, n) ~ H * H H — ^ . 

\(m + l) 2 (m + 2) 2 n 2 J 

Since Y^k=lW = TT> we concm de that (G^,dist c ) is a bounded metric 
space. 

Further, the resistance is bounded at infinity; or equivalently the voltage 
drop is "very" slow at infinity for the current flow induced by the experiment 
which inserts 1 amp at a particular place in G^ = No- 

The reason is that the conductance is "very" unbounded, or equivalently 
or more precisely, the resistance is O (n~ 2 ) for this particular (G, c) system. 
Some conclusions: The finite-energy solution v to (|5.60p is the function 
v : 7L 2 — > M, beginning with the values 0,-1/2, and —1 as follows: In 
Figs. 2-3 we list the values of v on the points in the interior square in 
{= Z 2 ) . The three values are prescribed in the centered square; and 
they then propagate into the quarter planes, with the value —1/2 in the 
NW and the SE quarter planes. 

6 Finite Dimensional Approximation 

6.1 Systems of Graphs 

Let G = (G^ ', G^) be an infinite graph satisfying the axioms from section 
2. In particular, we assume for every x in G^ that x itself is excluded from 
nbh(x); i.e., no x in can be connected to itself with a single edge. Let c 
any conductance function defined on and satisfying our usual axioms. 

In section [3] we showed that the corresponding Laplace operator A = A c 
is automatically essentially selfadjoint. By this we mean that when A is 
initially defined on the dense subspace T> (of all the real valued functions 
on G^ with finite support) in the Hilbert space H := £ 2 (G^), then the 
closure of the operator A is selfadjoint in TC, and so in particular it has a 
unique spectral resolution, determined by a projection valued measure on 
the Borel subsets the infinite half- line M+. 

In contrast, we note (Example IT. 1[) that the corresponding Laplace op- 
erator in the continuous case is not essential selfadjoint. 

This can be illustrated with A = —(d/dx) 2 on the domain T> of consist- 
ing of all C 2 -functions on the infinite half-line M+ which vanish with their 
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derivatives at the end points. Then the Hilbert space is L 2 (M+). 

So this is an instance where the analogy between the continuous case 
and the discrete case breaks down. 

In the study of infinite graphs G = (G^, CrW) and the corresponding 
Laplacians, it is useful to truncate and consider first a nested system of finite 
graphs G/v; then compute in the finite case and, in the end, take the limit as 
N — ► oo. Our approximation results here continue work started in |j0r77| , 
fj0r78l . 

Definition 6.1 In this section we prove specific results showing that the 
procedure works. While there are several candidates for designing the finite 
approximating graphs Gn = (G$ , G$ ) , we will concentrate here on the 
simplest: Starting with an infinite G = (G^°\G^), pick finite subsets of 
vertices as follows: 

Gf'cGfcGfc-cGPc-c^ (6.1) 

such that 

oo 

(Jg£)=g(°>. (6.2) 

N=l 

SetH : = £ 2 (G^), and H N = £ 2 (G^). Then the projection P N of H 
onto Tt onto TCn is multiplication by the indicator function XM°) > an d the 

projection onto the complement HQTIn is multiplication with X/M°hc where 

Ojv ) 

(G$) c = GW\G$ is the complement of G^ ] . 

The edges in Gn are simp 
lie in G^ ; i.e., if x,y G G^ , then 



The edges G$ in Gn are simple the edges in G, for which the vertices 



{xy) G G^ (xy) G G^ and x, y G G$ . (6.3) 

// a system (Gn)n<=n of graphs is given as in ip - (|fi.3|) . and if c : 
GW - R+ is a conductance function; we denote by cn the restriction of c 
to G$. 

Lemma 6.2 Let G = (G (0) , G (1) ) and c : -> R + be given as above. Let 
Gn be a system of graphs determined subject to conditions (|ff.ip - ([fi.gp . 
Let Atv be the graph Laplacian associated to (Gn,cn)- Then 

P N AP N = A N , for ViV G N. (6.4) 
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Proof. For v £ T> = finite linear combinations of {5 x \x £ G^}, we have 
(P N AP N v) (x) = xg n 0) c (xy) ((xg n v) (x) - (xg n v) (y)) 

y~x 

= ^2 cat (xy) (v (x) - v (y)) 

j/~rr in Gjv 

= (A N v) (x) ; 
proving the formula (|6.4p . ■ 

Lemma 6.3 Let G = (G^,G^), and c : G^ -> R + , be as in Lemma\EM 
and Definition \6.1\ Then for all v £ T> and x 6 G^ , we have the following 
formula for the difference operator A — An, N = 1, 2, ... : 

(Aw) (x) - (A N v) (x) = -XG% 0) c (^y) v (y) • ( 6 - 5 ) 

veG N 

In other words, the contribution to A — Ajy comes from the boundary of 
Gn = the edges e £ G^ s.t. one vertex in e is in Gpl and i/ie other in the 
complement. 

Proof. Using the previous lemma, we get 

(Av) (x) - (A N v) (x) = ^ (c (xy) - c N (xy)) (v (x) - v (y)) 

= -xg- n (x) c ( x y) v (y) • 

yeG N 

m 

Definition 6.4 Let G = (G<®,GM), and c : G^ R+ be given as in 
Theorem \ 5.1l and denote by A = A c the corresponding selfadjoint graph 
Laplacian. Setting 

POO 

S(t):= e- tx P(d\) (see (fOjl - fO^l ) (6.6) 
J o 

= e~ tA , <GR+; 

we see that 1 1 — > S (t) is a contractive semigroup of selfadjoint operators in 
£ 2 (G^); in particular, 

S(s + t) = S(s)S(t), Vs,t£R + and (6.7) 
S(0)=I P . 
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The semigroup consists of bounded operators while the infinitesimal gen- 
erator A = A c is typically unbounded, albeit with dense domain in £ 2 (G^). 
Moreover, the semigroup helps us identify dynamics as infinite graphs of 
resistors. 

Returning to approximations, as in Definition \6.1[ we now get a sequence 
of Laplacians An, N = 1,2, . . ., and a corresponding sequence of dynamical 
semigroups, Sn (t) ■ = e~* Ajv , N = 1,2, — 

Let N be fixed, and let 8Gm be the boundary of Gn (Definition \ 6.1\i . 
Then the finite matrix 

T N :={c{xy)) x y&dGN (6.8) 

is positive, and has a Perron- Frobenius eigenvalue Xn = Xn(PF) = the 
spectral radius o/Tjv- 

Theorem 6.5 Let (G, c) be a graph/ conductance system, and let (Gjv)iVeN 
ascending system of graphs such that (jfi.gp is satisfied. Let S (t) , and Sn (t), 
N = 1,2, . . ., be the corresponding semigroups of bounded operators. 
Then for all v G £ 2 (G^), we have the following estimate: 

\\S(t) v-S N (t) v\\ £2 < X N (PF) t \\v\\ £2 , Vt G R + , N = 1,2, ... . (6.9) 

Proof. With the use of ( \5.8\) - ( 15.9)) , we get the integral formula: 

e ~tA N _ e -tA = f* e -(ts)A (A _ A ^ e -sA N dg ^ 1Q) 

Jo 

Since the operators on both sides in (|6.10p are bounded, it is enough to 
verify the estimate (|6.9p for vectors v in the dense domain D. 

Using new Lemma 16. 3) we get the following estimates on the respective 
f 2 -norms: 

\\S(t)v-S N (t)v\\ e2 < f || (A - A N )S N (s)v\\ £2 (bydHUD) 

J o 

<X N {PF) / ||5jv(s)u|| ds (by LemmaESand flgSJ)) 
Jo 

< X N (PF) \\v\\ t 2 [ ds 
Jo 

= X N {PF)t\\v\\ £2 , 
which is the desired conclusion. ■ 
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6.2 Periodic boundary conditions 



Example 6.6 We now compare Example \5.i$ with an associated family of 
finite graphs G N where N G N. Let Z N = Z/NZ ~ {0, 1, 2, . . . , N - 1} be 
the cyclic group of order N. Introduce nearest neighbors as in Example \5.2\ 
(the TL-case) with the modification for Gn given by ~ (N — 1), in other 
words that there is an edge connecting to N — 1. 



It follows that the graph Laplacian Ajv for Gn is the given by the finite 
matrix 






-1 













2 -1 
-1 2 
-1 



-1 





-1 

2 



The spectrum of Ajy is as follows: 

'2vr/c\ 



spec (Ajv) 



2 1 



cos 



irk 



4 sin' | - 



N) 
\k = 0,1, 



0,1,... ,JV 
,N-1 



(6.11) 



Comparing with (j5.34j) - (|5.35[) . we see that the spectra converge in a natural 
sense; with the infinite model in Ex l5.2l being a limit of iV-periodic boundary 
condition as N — > oo. 
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Fig. 4. The vertices and edges in Gn for N = 3, 4 and 5. 

The spectrum of the cyclic graph Laplacian A at of the graphs Gn, N = 
3,4 and 6 (in Fig. 3) may have multiplicity; see (16. lip . This holds in fact 
for all values of N. Specifically, 

spec (A3) = {0, 3} with A = 3 having multiplicity 2. 

spec (A4) = {0, 2, 4} with A = 2 having multiplicity 2. 

spec(Ag) = {0, 1,3,4} now with A = 1 and A = 3 each having multi- 
plicity 2. 

Hence for D = 1, we get the following distinction between the spectral 
theory of the cyclic graph Laplacians Atv for N < 00 on the one hand and 
A in Example 15.21 on the other: The commutant of A is an abelian algebra 
of operators in £ 2 (Z), while the operators in £ 2 (Zjv) which commute with 
Atv form a non-abelian algebra. 

Proposition 6.7 (Cyclic graphs) Lei N G N ; N > 3; and let Gn be the 

corresponding cyclic graph with graph Laplacian An; i.e., with 

G^ = {0,1,2,. ..,N-l}. 
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Then the voltage potential v G £ 2 (Z/v) solving Ajyv = 5q — 5\is 

( v = 




< 



VN-2 = —j? 

VN-l = -jv- 



Proof. A direct computation; see also Fig. 4, and eq. (|6.10p . ■ 

7 Boundary Conditions 

In the study of infinite graphs G, boundary conditions play an important 
role; for example if a current escapes to infinity in "finite time," conditions 
must then be assigned "at infinity." 

One way to do this is to first do computations in a system of finite graphs 
G n which exhausts the given graph G in a suitable way. Do computations on 
each finite subgraph Gn of the fixed infinite graph G, and then take the limit 
as N tends to infinity. There are several ways one may do the computations 
on each individual Gn , for example look for symmetry, or look for a suitable 
periodicity, or similarity up to scale. In the simplest cases, this allows the 
use of a finite Fourier transform, thus making Gm periodic, or cyclic. The 
case of G = Z (the rank-1 integer graph), and Gn = the cyclic group of 
order TV" is done in all detail in Example 16.61 above. 

Some advantages of the cyclic approach: One, the spectrum comes out 
given explicitly by a closed formula, thus making it clear how the limit 
TV" — > oo works also for spectra, getting the continuous spectrum in the 
infinite limit. 

Example 7.1 In this section we compare the two cases, continuous vs. dis- 
crete. As noted, our graph Laplacians are second order (or more than second 
order) difference operators in a generalized sense. 

They have spectrum contained in the half-line [0,oo), so generalizing 



V:= {veC 2 (0, oo) \v, v', v" G L 2 (0, oo) , and v (0) = v (0) = 0} ; (7.2) 




(7.1) 



with the Hilbert space TL : = L 2 (0, oo), and dense domain 
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i.e., with vanishing boundary conditions on v and v 1 (x) = ^ at x = 0. 
We get the spectral estimate: 

(v, Av) L 2 > 0, \/v £ V. (7.3) 

A simple verification shows that for the adjoint operator A* we have: 

dom (A*) = {veL 2 (0, oo) \v', v" G L 2 (0, oo)} . (7.4) 

Comparing (|7.2p and (|7.4p we see that A results from A* by "removing" 
the two boundary conditions which specify the domain T> of A. 
Moreover, the defect space 

V + : = {v G dom (A*) |A*u = -v} (7.5) 

is one-dimensional; in fact, 

V + = Ce~ x . (7.6) 

The selfadjoint extensions of A on T> are parametrized by pairs of num- 
bers A, B G R, not both zero, such that 

Av (0) + Bv' (0) = 0. (7.7) 

Example 7.2 Let G = (G^ q \G^) be the following graph generalizing the 
continuous example: 

: = No = {0,1,2,3,...}, (7.8) 
:={(01),(n,n±l); n G N} . 

Pick A > 1, and set 

c(n,n + 1) = A n+1 . (7.9) 
Then the corresponding graph Laplacian is unbounded; and 

(Av) (0) = Xv - Xv i; (7.10) 

(Av) (n) = -X n v n . 1 + A" (1 + A) v n - A n+1 t» n+1; Vn G N. (7.11) 
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For domain T>, we take all v G £ 2 (No) s.t. v n = except for a finite set 
of values of n. the matrix representation of A is presented in Fig. 4 : 

( A -A ••• ••• \ 

-A A(l + A) -A 2 

-A 2 A 2 (l + A) -A 3 




vn+l 



; -A n 

O -A n A n (l + A) 

\ -A n+1 , 

Fig. I 

By Parseval's formula, we have the isometric isomorphism £ 2 (No) 
7i + = the Hardy space of analytic functions on D = {z E C; \z\ < 1} 



oo 

r !::):= ^ /'„ . 



n=0 

and 



m\h + = El^l 2 - ( 7 - 12 ) 

n=0 



7n the Hardy space representation we have 

(At;) (z) = (1 + A) « (Az) - Xzv (Xz) - z~ l v {Xz) (7.13) 

on the dense space of functions v on C which extend analytically to D\ : = 
{z G C; |z| < A}. 

We now show that there are no non-zero solutions to 

A* x v = -v, (7.14) 

i.e., v G dom(A^); equivalently V + (A) = {0}; the defect space for the 
operator A\ is trivial. So this is a direct verification that A\ is essentially 
selfadjoint; and contrasting with ( |7. 6\) above. 

To see this, combine (|7.13p and (|7.14p . It follows that every solution v 
to (17.141) must have an infinite-product representation given by 

. . (z — 1) (Xz — 1) . . . H 

v[z)= y - ^ L v(Xz); 7.15 

Az 
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and the limit of finite products as follows 

n—l „ 

{z- 1) n (A fc 2- 1) (X n z- 1) 

k=l 



These products do not have a non-zero representation consistent with the 
isomorphism (|7.12j) . and with (|7.12j) . 



8 Appendix 

A Heisenberg's Infinite Banded Matrices 

We proved in sections 3 through 5 that in general, graph Laplacians A c 
are essentially selfadjoint operators in the £ 2 sequence-Hilbert space. Recall 
that the axioms for our graph Laplacians include the following given data: 
A graph G = (G^°\ G^) and a fixed positive conductance function c defined 
on the set of edges G^ 1 '. Every vertex x of G is connected to a finite set 
of neighbors in G^ . For every fixed x in G^ , this implies finiteness of the 
set of y in G^ for which c(xy) is nonzero. This means in turn that the 
natural matrix representation of the operator A c is banded; see section [5] for 
the Definition. Note however that we place no boundedness restrictions on 
the conductance function c. 

Our proof of essentially selfadjoint for the operator A c uses this hand- 
edness property in an essential way. In fact, starting with an infinite by 
infinite matrix, it is generally difficult to turn it into a linear operator in a 
Hilbert space unless it is assumed banded, see section [H and the references 
cited there. 

The purpose of this section is three-fold. 

First to make precise the operator theory of banded infinite by infinite 
matrices; and second to show that the infinite matrices used in represent- 
ing the operator algebra generated by Heisenberg's quantum mechanical 
momentum and position observables consists of (infinite) banded matrices. 
Thirdly, we use Heisenberg's (and Born's) computations to exhibit such 
banded operators which are not essentially selfadjoint. The simplest such 
matrix M is as follows: let P be Heisenberg's momentum operator and Q 
the (dual) position operator. Then we show that the monomial M = QPQ 
is banded, but not essentially selfadjoint. In fact, its deficiency indices are 
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Definition A.l Let L be a countable (typically infinite) set, and let m : 

L x L — * C be a function on L x L. We say that m is banded iff for every 
x G L, the set 

{y G L\m(x,y) / 0} (A.l) 

is finite. 

Let £ 2 (L) be the sequence space with norm 

\\v\\% :=^2\v(x)\ 2 <oo. (A.2) 

The sum on the right is the supremum of all the numbers X^eF \ v ( x )\ 2 
as F ranges over all finite subsets in L. 

Let T> be the dense subspace of all functions v : L — > C such that the 
support set 

{x G L\v (x) / 0} (A.3) 
is finite. Equivalently, setting 

the space D is then the linear span of the set of functions {5 X \x £ L}; 
and these functions form an orthonormal basis for £ 2 (L). Moreover, every 
Hilbert space TL is isomorphic to £ 2 (L) for some set L. The set L is countable 
if and only if 7i is separable. 

Lemma A.2 Let m : L x L — > C be a banded function. For v G T> C £ 2 (L), 

(Mv) (x) = m (x, y) u (y) . (A. 5) 

yeL 

Then M defines a linear operator M : T> — > T>, with a well defined adjoint 
operator M* . Moreover, 

DCdom(Af) (A.6) 
where dom(M*) is the domain of M* . 

Proof. When x G L is fixed, the sum in (|A.5P is finite because the set (jA.ip 
is finite by assumption. Using finiteness of both sets (|A.1|) and (IA.3j) we 
conclude that Mu in (|A.5h is in P if w is. And so, in particular, Mv G ^ 2 (L); 
see (|X2l) and (hOl) . 
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To establish the inclusion "C" in (|A.6[) . we must show that for every 
v £ T>, there is a constant K = K (v ) such that the following estimate 
holds: 

\{Mu,v)p\<K\\u\\ p , forV?ieP. (A.7) 
The expression on the left in (|A.7j) is 

^2^2<m(x,y)u(y)v(x) . (A.8) 

x,y£L 

But the terms in this double-sum vanish outside a finite subset in L x L an 
account of assumptions (jA.lj) and (IA.3j) . 

The modulus-square of the sum in (|A.8P is estimated by Schwarz by: 



y£L yeL 



^ j m(x,y)v(x) 



which yields the desired estimate (|A.7p . ■ 

Corollary A. 3 Let M be a linear operator in a Hilbert space H. Then M 
has a banded matrix representation if and only if there is an orthonormal 
basis (ONB) in Tt, {e x \x £ L] such that the linear space V spanned by 
(ex)xeL * s mapped into itself by M . 

Corollary A. 4 In that case the matrix entries of M are indexed by L x L 
as follows: 

m(x,y) := (e x ,Me y ). (A.9) 

Proof. Only the conclusion (|A.9|) is not contained in the lemma. Now 
suppose some operator M in 7i satisfies the conditions, and let (e x ) xeL be the 
associated ONB. Then Me y G H ~ I 2 (L), so Me y = ^{e^Me^e^, 
and 

\\Me y \\ 2 n = ^2\(e x ,Me y }\ 2 (A.10) 

holds by Parseval's formula. The conclusion (|A.9h follows. ■ 
Corollary A.5 Let G = (G^,G^) and 

c : G (1) -> R + 



47 



be a graph system satisfying the axioms in section^ Let {S x \x G G^} be the 
canonical ONB in £ 2 (G^). Then the graph Laplacian has a corresponding 
banded matrix representation as follows: 



(5 x ,A c 6 y ) = 
Proof. Recall the function 



— c (xy) if y / x and y ~ x 

B c 0) ify = x 

if ' y <]L x and y ^ x. 



B c {x) 

on the right-hand side in (jA.lip . 
Since, for v G V, we have 



(A.ll) 



(A.12) 



(A c v) (x) : = (xy) (t> (x) - v (y)) 
setting v = 5 y , we get 

(A c <5 y ) (x) 



(A.13) 



B c (x) if y = x 
— c (xy) if y ~ x 
if y </■ x 



from which the desired formula (|A,lip follows. ■ 

Heisenberg introduced oo x oo matrix representations for the operators 
of momentum P and position Q in quantum mechanics. 

In the simplest case of one degree of freedom, they are as follows: 



/ 1 

1 V2 

V2 

\/3 




































n 









Vn+l 









V : 
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and 



1 

2i 



( -1 

1 -V2 ••• 

y/2 -V^2 

: : : -y 7 ^ 7 ! 

yfn~^l 

y/n + T 



n 





Set N : = {0, 1,2, . . .} = Z+ U {0}, and W : = ^ 2 (N ). Then the two 
matrices P and Q are represented by the following second order difference 
operators, having the same form as our graph Laplacians (1A.13P - 

1 



(Pv) (n) = - (Vn - 1 v(n-l) + y/nv(n + l)); (A. 14) 

and 

(Q V ) ( n ) = — (v 7 ^ 7 ! w (n - 1) - y/n v (n+ 1)) , (A.15) 
2i 

for Vi* 6 P, Vn £ No; where i = 

It is well known that both P and Q, as in (|A.14p and (|A. 15[) . are essen- 
tially selfadjoint. 

It follows by the above lemma that 

M: = QPQ (A.16) 

is also a banded operator., referring to the canonical ONB {e n |n 6 No} in 
^ 2 (N ). 

Caution: All the operators P, Q, and M are unbounded, but densely 
defined; see |j0r77| , |J0r78| , |Sto51| . 



Proposition A. 6 The operator M in <\A. 16\i is Hermitian, and has defi- 
ciency indices (1, 1); in particular is not essentially selfadjoint. In fact, it 
has many selfadjoint extensions; a one-parameter family indexed by T. 

Proof. By the Stone-von Neumann uniqueness theorem, the two operators 
P and Q in (|A.14p and (|A.15P are unitarily equivalent to the following pair 
in the Hilbert space I? (R) of all square-integrable functions on the red line: 

(Pf)(x) = -^-f(x), (A.17) 

i dx 
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and 

(Qf) (z) = xf (x) , for V/ G L 2 (R) , x G R. (A.18) 

For domain £> in (|A7T7|) and ([A~T8j) . we may take T> : = (R) , or the span 
of the Hermite functions. 

From the representations (|A.14p - (|A.15|) . it follows that the operator M : 
= QPQ in (|A.16p commutes with a conjugation in the Hilbert space; and 
so by von Neumann's theorem (see Remark 14. 2p . it has deficiency indices 
(n, n). We will show that n = 1. Hence we must show that each of the 
equations M*v± = ±iv± has a one-dimensional solution space in H. 

Taking advantage of Schrodinger's representation ()A.17P - (|A.18p . we ar- 
rive at the corresponding pair of ODEs in L 2 (R): 

x-^(xf) = ±f(x). (A.19) 

By symmetry, we need only to treat the first one. 
A direct integration shows that 

f(x) = < ltx > (A.20) 

V 7 \ if x < V ' 

solves (IA.19P in the case of "+" on the right hand side. Also note that 
(|A.20p is meaningful as all the derivatives of x _1 exp(— -) for x G R+ tend 
to when x — > 0+. This means that the two separate expressions on the 
right-hand side in ()A.20p "patch" together differently at x = 0. 

By the reasoning alone, we conclude that M has indices (1,1). As a result 
of von Neumann's extension theory, the distinct selfadjoint extensions of M 
are then indexed by T = {z G C| \z\ = 1}. If z G T, and if f± are normalized 
solutions to (|A.19p . then the extension M z is determined by 

M t (J + + zf-)=i(U-zf.). 



Example A. 7 Let P and Q be the canonical momentum and position op- 
erators; see ( \A.14\) - f\A.15§ , and let 

H : = P 2 - Q 4 (A.21) 

be the Hamiltonian of a "particle-wave" in one degree of freedom, corre- 
sponding to a repulsive x 4 potential. Then the reasoning from above shows 
that H is a banded ooxoo matrix. As an operator in I 2 (Z), H has deficiency 
indices (2, 2). 
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